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PREFACE 

This Is Volume I of a two volume final report for NASA 
Contract No. NAS8-32171 entitled, "SRB Reentry Acoustic Environ- 
ment and Flexible Heat Shield Similarity Verification." In this 
volume, we report on those aspects of the program bearing upon 
the aeroelastic stability of the SRB heat shield during reentry. 

Our assessment of the reentry acoustic environment, along with a 
derivation and discussion of the relevant scaling parameters, is 
presented in Volume II. 

The authors would like to acknowledge many helpful discussions 
with W.W. Clever and E.H. Walker of NASA MSFC. Their cooperation 
and assistance throughout the program is appreciated. 


PAOB BUNK NOT FlTJfl® 
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I. INTRODUCTION 

The heat shield on the Solid Rocket Booster (SRB) Is a fabric 
curtain that covers the annular opening between the lower edge of 
the aft skirt and the glmbeled rocket nozzle. Its primary purpose 
is to protect components within the cavity between the skirt and 
the nozzle from heat radiated by the exhaust plume during ascent, 
but it also reduces fluctuating pressure levels within the same 
cavity during reentry. A schematic drawing of the aft end of the 
SRB including the heat shield is shown in Pig. 1. 

The SRB tumbles after being separated from the Orbiter, and 
then settles into a tail-down attitude as it reenters the atmos- 
phere. Maximum dynamic pressure is reached when the reentry 
Mach number is about 2.75. During this critical phase of reentry, 
the angle between the longitudinal axis of the booster and the 
relative wind varies between 110° and 180°, with the aft end 
falling foremost. The heat shield, is thus directly exposed to a 
suoersonic flow at high dynamic pressure. 

Wind tunnel tests of a 3? scale model of the aft portion of 
the SRB equipped with partially scaled heat shields were conducted 
for the purpose of measuring fluctuating pressure levels in the 
aft skirt region.* During these tests, the heat shields were 
observed to oscillate violently, the oscillations in some instance 
causing the heat shields to fail. High speed films taken during 
the tests reveal a regular pattern of waves in the fabric starting 


The mass/unit area of the model heat shields bore the proper 
relationship to full scale. Other similarity parameters bearing 
upon forced or self-excited aeroelastic response were not scaled. 
Mass-ratio scaling is appropriate for modelling the "mass law" 
response of an essentially limp membrane to enforced fluctuating 
pressure loads. 
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near the flow stagnation point and progressing around both sides 
of the annulus. The amplitude of the waves was too great, and 
their pattern too regular, for them to be attributed to the 
fluctuating pressure levels measured during the tests. 

The purpose of the investigation reported here is to deter- 
mine the cause of the oscillations observed in the model heat 
shields, and to determine whether or not similar oscillations 
will occur in the full scale SRB heat shield during reentry. 
Suggestions for modifying the heat shield so as to avoid the 
oscillations are provided, and recommendations are made for a 
program of vibration and wind tunnel tests of reduced-scale 
aeroelastic models of the neat shield. The recommended program 
provides an enhanced level of confidence in establishing a heat 
shield design that will survive reentry and provide protection 
for components contained within the aft skirt cavity. 


3 
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2. QUASI-STEAOY AERODYNAMIC ENVIRONMENT DURING REENTRY 

Data on the expected reentry trajectories of the SRB have 
been provlded%y NASA Marshall Space Plight Center (Gallaboff, 
1976). Of most Importance to the present investigation are the 
piedlcted maximum values of dynamic pressure and the predicted 
attitude or angle of attack of the SRB during reentry. The 
maximum dynamic pressure (to be exceeded during only 51 of re- 
entries) is shown In Pig. 2 as a function of reentry Mach num- 
ber M^. The greatest dynamic pressure is encountered when the 
reentry Mach number Is about 2.75. 

During the portion of reent: •' at which maximum dynamic pres- 
sure is encountered, each SRB is expected to be executing a con- 
ing motion with an angle of attack greater than 90°. That Is, 
the boosters fall nozzle-end first (see sketch in Fig. 3). The 
angle of attack Is not constant, but varies randomly about an 
equilibrium value. The angle of attack is expected to fall out- 
side the range shown in Fig. 3 during only 5/K of all SRB re- 
entries. 

Even though the angle of attack Is basically a random vari- 
able during reentry, the rate of variation of angle of atta k 
is very small compared to the frequencies of the acoustic and 
aeroelastic disturbances being addressed in this report, so the 
flow around the aft end of the SRB can be assumed to be quasi- 
steady; that is, identical to the steady flow that would exist 
at the same instantaneous angle of attack. We can therefore 
infer the basic features of the flow over the aft skirt region 
of the SRB from th 5 Schlleren photographs and time-averaged 
pressure data recorded during wind tunnel testing of a 3 % 
scale model in the 1.02m x 1.02ra (40 in x 40 in) Supersonic Wind 
Tunnel (A) at the Arnold Engineering Development Center (AEDC). 


4 











BBN Report No. 3532 


Bolt Beranek and Newman Inc. 



FIG. 2. DYNAMIC PRESSURE VS MACH NUMBER 
DURING REENTRY 










Report No. 3532 Bolt Beranek and Newnan Inc. 


These tests were conducted at Mach numbers of 2.75 and 3.5. 
The angle of attack range was 110° to 139** • The test program 
was designed to d3termlne maximum acoustic loads, so the test 
procedure was to record only rms pressure levels at selected 
points In the aft skirt region over the entire angle of attack 
range. Schlleren photographs were taken only at thos^ angles 
of attack which produced the highest rms pressure i luctuatlons, 
so the shock structure Is not known at all angles of attack. 
However, the available Schlleren of photographs show the situa- 
tion depicted In Fig. 4. 

At the angles of attack at which the greatest aeroelastlc 
response and acoustic loads occurred an oblique shock wave was 
attached to the lip of the aft skirt at the stagnation point. 

In all cases checked, the angle formed by this shock wave with 
the direction of flow agreed (within ±1°) with that predicted by 
Invlscld compressible flow theory for the two-dimensional flow 
over a wedge having the flow deflection angle 0 shown In the 
figure. This was Interpreted as evidence that, at least In the 
plane of lateral symmetry, the flow was essentially a 2-D wedge 
flow. To confirm this, time-averaged (DC) pressures recorded 
at several locations on the face of the aft skirt were plotted 
vs. angle of attack. Typical results are shown In Pigs. 5 and 
6 (M^ » 2.75 and 3.5, respectively). In each figure, the solid 
line Is the pressure that would exist on the face of a wedge 
that turned the flow through the angle 0 * a - 90°. The lines 
are extended only to the angle at which the shock would detach 
from the leading edge In a 2-D flow. At higher angles of attack, 
the shock will stand off from the skirt, and the flow will be 
subsonic over some or all of the face of the skirt. 
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The locations of transducers 1 through are shown in Pig. 7* 
(locations PI through P<1). Transducer #2 is located immediately 
upftream of a small lip extending around the periphery of the 
tTancated rocket nozzle. At high angles of attack, this trans- 
ducer would be expected to read stagnation pressure, as indeed 
it is apparently doing in both figures for o > 120®. The other 
transducers, however, show reasonable agreement with the wedge 
f 1 *w pressure curve for all angles less than the shock detach- 
mc it angle. 

In the following sections, we will assume that mean flow 
conditions over the face of the heat shield and the skirt proper 
can be related to free stream conditions by computing the flow 
over a wedge at the appropriate upstream Mach number and turning 
angle. The requisite theory and formulae are given in Liepman 
and Roshko (1957) and in NACA Report 1135 (1953)* and will not be 
repeated here. For the present, we show in Figs. 8 and 9 the 
local Mach number and local dynamic pressure (normalized by the 
'ree stream value) as functions of angle of attack and free 
stream Mach number. These figures apply equally to any geo- 
metrically scaled wind tunnel model, and to the SRB Itself. 

For the model, of course, M and q denote wind tunnel test 
* ^00 *00 

section conditions, whereas at full scale they are the reentry 
Mach number and dynamic pressure. 

To unde”s uand the aeroelastlc oscillations of the heat 
shield, it is important to understand the following two points, 
implicit in Figs. 8 and 9: 


♦Figure taken from Porter (1976). 
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(1) The local Mach number over the heat shield Is always 
less than or equal to the free stream Mach number. 

As the angle of attack Increases from 90®, the local 
Mach number drops from the free stream value to 1.0. 
At still higher angles of attack, the flow over some 
portion or all of the heat shield is subsonic. 

(2) The local dynamic pressure may exceed the free stream 
value, in spite of the viscous losses that occur in 
the oblique shock. By comparing the two figures, it 
can be seen that even where M^ « 3.5, the greatest 
local dynamic pressures occur at a local Mach number 
between 1.0 and 2.0. 

Anticipating later results, the aeroelastic response of 
the heat shield is expected to be moat severe w ,n the local 
Mach number exceeds 1. As shown in Fig. 10, the local Mach 
number is expected to exceed 1 only when the reentry Mach 
number is greater than 2.6. This potentially troublesome 
regime could he eliminated entirely if the expected angle 
of attack range could he narrowed for descent Mach numbers 
greater than 2.6. To eliminate locally sypersonic velocities 
over the heat shield during the entire reentry, it would be 
necessary to insure that the angle of attack exceeds 132® at 
M^ = 2.75 and 137° at M^ = 3*5. Since the local Mach number 



is very sensitive to changes in angle of attack near 

the angles quoted above Incorporate a safety factor of 10®. 


is, the theory predicts that = 1 when a = 122® for 

M * 2.75, and when a * 127° for M = 3.5* 


= 1 , 
That 


Elimination of locally supersonic flow over the heat shield 
will not preclude the possibility of unacceptable heat shield 
oscillations, but it will reduce their severity. The addltlcnal 
structural stiffness required to eliminate the oscillations will 
be correspondingly reduced. 
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3. MODEL HEAT SHIELD OSCILLATIONS OBSERVED IN WIND TUNNEL TESTS 

The characteristics of the heat shield oscillations observed 
during the AEDC tests are summarized In Table 1. Anticipating 
later developments, the corresponding results from flutter theory 
are also shown. The entries in this table indicate the ranges 
encountered for the quantities listed; no attempt has been made 
to distinguish between specific heat shield designs* or flow 
conditions (Mach number, angle of attack, tunnel density, etc.). 
No instrumentation intended specifically to measure the heat 
motion was installed on the model, so only the high speed film 
provides any specific information about the oscillations. These 
high speed films were not taken on all runs, but all films which 
display the oscillations shov; a regular wave motion starting at 
the portion of the annular heat shield nearest the stagnation 
point, and progressing around both sides of the annulus in the 
direction of air flow. The wave motion is symmetric about the 
diameter containing the stagnation point, with pairs of waves 
moving around the annulus from the stagnation point and meeting 
one another 180° from their point of origin. No reflected waves 
were discernible. 

TABLE 1. MODEL HEAT OSCILLATIONS 



Measured 

Flutter Theory 

Wavespeed 

152-178 cm/sec 

150-760 cm/sec 

Wavelength 

4-12 cm 

7-37 cm 

Frequency 

120-250 Hz 

50-800 Hz 

Amplitude 

1.2 cm 

( approximately ) 

— 


*Heat shields of three different materials and area mass densities 
were tested. 
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No forced-response phenomena known ro the authors can explain 
this behavior. Even If a source of excitation at a single 
frequency were present In the flow, traveling waves would not 
be indicated unless tne disturbance was limited to one region 
of the heat shield (specifically, the region from which the 
waves originated), or the disturbance Itself were convected 
around the annulus with the observed speed of the waves in the 
heat shield. Note that the speed of the waves is consld'^rably 
less than the fluid flow velocity which Is on the order of 2500 
cm/sec. Thus, the waves are not likely to be a passive response 
to a vortex or other disturbance convected with the flow. 

On the other hand, the flutter theory presented in Sec. ^ 
predicts the following features of the motion observed In the 
wind tunnel tests: 

(1) The basic wave-like character of the oscillations; 

(2) Frequencies and wavelengths which are In fair quanti- 
tative agreement with experiment; and 

(3) Good agreement with experiment as to the range of 
angles of attack at which the aeroelastic response 
Is most severe. 

Two additional points of agreement should be mentioned here. 
Although no attempt has been made to predict the amplitude of 
the heat shield oscillation, nonlinear panel flutter theory pre- 
dicts that it should be on the order of the sag radius (Dowell, 
1970). The sag radius of the model heat shields was R = 1.1^4 cm 
(0.il5 in), which Is roughly the amplitude of the oscillations 
seen In the films. Nonlinear flutter theory also predicts that 
a static pressure differential should reduce the amplitude of 
the response when the pressure is higher on the concave side of 
the heat shield, and Increase It when the pressure Is higher 


18 




I I 1 I I I „J 


Report No. 3532 Bolt Beranek and Newma 


on the convex side (Dowell, 1969 and 1970). Detailed calcula- | 

tlons of the effect of a pressure differential are discussed | 

In Sec. here we mention only that one test run In which the | 

pressure within the cavity behind the heat shield was greater | 

than that outside produced violent oscillations and Immediate ; 

failure of the heat shield. 1 

*1 

The success of the flutter theory In predicting the ^ oIc 
features of the heat shield oscillations and the flow c; Ions 

under which they occur lead us to conclude that the osci ’ons 

were caused by aeroelastlc flutter, rather than some forced 
response phenomenon. 



n Inc. 




N. 
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4. HEAT SHIELD FLUTTER: BASIC THEORY AND APPLICATION TO 3* 

SCALE MODELS 

In this section, we present basic theoretical results for 
heat shield flutter and discuss their application to the 3 ? scale 
heat shields used in the AEDC tests. Two separate theories are 
presented, ore for low supersonic speeds and another for high 
supersonic speeds, because of basic differences in the parametric 
dependence of the a-'roelastlc response in these two flow regimes. 
The funJament£ ’ difference between the two theories Ir the method 
by which the aerodynamic loads on the oscillating heat shield are 
calculated. The simplified piston theory expression used in the 
high supersonic theory is usually considered to be acceptable 
for local Mach numbers greater than 1.5 or 2.0 (see Dixon ri066] 
or Dowell [1972]). In this report, we have chosen M t 1.5 as 
the division between the two regimes. Further discussion of 
this point is provided in Sec. 4.3. 

4.1 Prediction of Aeroelastic Loads 

Under most test conditions, the oscillations encountered dur- 
ing the AEDC tests did not cause Immediate failure of the heat 
shields. Similar oscillations in the full scale heat shield might 
also be acceptable, if it could be ascertained that they would not 
cause premature failure of the heat shield or nearby components 
during reentry. This approach would require that stresses and an- 
chorage loads caused by flutter be determined with an accept- 
able level of confidence. The relevant theory is inherently 
nonlinear (Dowell, 1970), and requires lengthy numerical com- 
putations to treat even simple structural configurations (such 
as, for example, a rectangular homogeneous elastic plate with 
no applied structural loads). A r*»cent a tempt to compute 
etveeBee in a fluttering panel exposed to In-plane loads 
and a static pressure differential met with limited success 
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(Ventres and Kang, 1973)* In our opinion, predicting aeroelastlc 
loads by theory alon* Is beyond the state of the art. Therefore, 
we have limited our attention to predlCv.ng the onset of flutter, 
with the ultimate objective of using theory and experiment* to 
design a heat shield that will not flutter during reentry. 

4.2 Low Supersonic Theory 

The hlg* speed films taken during the AEDC tests show waves 
progressing around the circumference of the annular heat shield. 
The wavelength of these waves Is comparable to the width of the 
annulus (more precisely, the difference between the inner and 
outer radii). Since the wavelength Is much smaller than either 
the semi-circumference or the mean radius of the heat shield, 
we have simplified the geometry for analytical purposes by 
"straightening out" the heat shield. That Is, the aeroelastlc 
model is not an annulus but an infinitely long panel having the 
same width as the a..nulus. The airflow is assumed to be parallel 
to th<- side edges of this infinite panel, as shown in Fig. 11. 

To justify this flutter model, a series of companion cal- 
culations were made for a set of aluminum panels of various geo- 
metries. The quantity calculated In all cases is the thickness 
required of a homogeneous alumimun plate having the same width 
as the full scale SRB heat shield, so as to prevent flutter 
during reentry. The results are summarized In Table 2. They 
demonstrate the following points: 

(1) The thickness required to prevent flutter Is greater 
when the flow is circumferential. 


*Recommended tests are discussed in Sec. 6. 
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TABLE 2. THICKNESS REQUIRED TO PREVENT FLUTTER OF A FULL SCALE 

HEAT SHIELD C0NSISTIN6 OF A FLAT ALUMINUM PLATE 76.2 cm 
(30 In) HIDE: TEST OF ASSUMPTIONS ON CRITICAL DIRECTION 

OF FLOW AND LENGTH OF EQUIVALENT STRAIGHTENED HEAT SHIELD. 
LOCAL MACH NUMBER -1.3. 


Direction of Flow 



Lengthwise 

Across Hidth 

Plate Length 

(Circumferential ) 

(Radial } 

76.2 cm 

0.79 cm* 


(30 in) 




0.86 cm** 

0.72 cm*** 


* Prom Pig. 2 of Dixon (1966). The mean semi-circumference of 
the full scale heat shield is 710 cm. 

** Using present low supersonic traveling wave ^heory based on 
Dowell (1966). 

*»* Prom Dowell (1975). Assumes * 25- 


23 





T 


% 

i 


1 


1 


Report No. 3532 


Bolt Beranek and Newnan Inc. 


ti! 


(2) The effect of finite length in the direction of air- 
flow is small, and in any event, assuming an in- 
finitely long panel is conservative (i.e., requires 
a thicker plate to prevent flutter). 

A schematic drawing of the annular heat shield and the 
straightened out approximation are shown in Pig. 11. The full 
scale heat shield is constructed so that the heat shield fabric 
forms a sag radius of 37.3 cm (1^1.7 in) between the nozzle 
anchorage ring and the aft skirt. The equivalent radius (at 
3X scale) was formed into the model heat shields, and is shown 
in Pig. 11 as the sag radius R. This has been incorporated in 
the straightened model as cross-stream curvature of the same 
radius R. 

The analysis that follows is basically an adaptation of the 
method devised by E.H. Dowell (1966) to treat the flutter of an 
infinitely long flat plate. His analysis has been extended to 
Include cross-stream curvature and a static pressure differential 
(unequal pressures on the inner and outer surfaces of the heat 
shield). 

The linearized equations of motion of a plate curved in 
one direction are: 

DV**w + I ^ - N — - N 

” dx* "" 3x* y 3y* 


+ m D 

3t» 


7>-f ^ 

^ 3x* 


J 


.1 


1 
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As depicted in Fig. 11, the x-axls extends along the Infinite 
dimension of the plate, whose side edges are located at y > ±b/2. 
The airflow Is In the direction of Increasing x. Refer to the 
List of Symbols for other definitions. 

In Eq. (1), and are the middle surface rtress resultants 
caused by the static pressure differential Ap. Since the sag 
radius R Is roughly equal to b/2, the heat shield takes on the 
shape of a half-torus (bisected at the meridian plane). There- 
fore, '.he stresses In the heat shield due to the pressure dif- 
ferential will be the same as those in an Internally pressurized 
torus (or Inner tube), which are approximately* 

n . 


N * ApR 

y 


The method of solution is to assume a traveling wave solu- 
tion for the normal deflection w(x,y,t) and the Airy Stress 
function iKxjyjt):** 


io(ct-x) 


w = Re 


ia(ct-x) . 
<» « Re < e > 


if) 


*Exact formulas are provided In Dyra (197^). The expressions given 
above become exact when the ratio of the sag radius R to the 
circumference of the heat shield vanishes. 

**These expressions Imply sln?)ly supported boundary conditions 
on w and stress-free boundary conditions on i|». 
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These expressions are substituted Into Eq. (2) to determine 
♦ in terms of w; 


♦ 



(5) 


If we now substitute Eq. (4) into Eq. (1), multiply by cos(iry/b) 
and integrate from y ■ -b/2 to y » +b/2, and finally use Eq. (5) 
to eliminate we obtain 


V b^ / r2 




N o* + N — 
X y ^2 


2 2 

- me 


b '' 

2 w e 


ia(ct-x) 


b/2 


j p cos 1^) dy = 0 


( 6 ) 


-b/2 


The last term in Eq. (6) contains an Integral of the pressure p 
caused by the travelling wave deflection w. In his investigation 
of the flutter of Infinitely long flat plates, Dowell (1966) re- 
duced this expression to a function F(n) of a single variable n, 
which he evaluated numerically: 

+b/2 

I p cos(iry/b)dy * -pa* 

-b/2 


(U - c)* ^ F(ri) 

ia(ct-x) 


where 
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The function F(n) is defined In Dowell (1966) and Is displayed 
graphically In Fig. 2 of Dowell (196^1). Equation (7) has been 
used for the aerodynamic pressure on the heat shield. This 
amounts to assuming that the relation between the heat shield 
motion and the resulting pressure is not substantially changed 
by the cross-stream curvature (heat shield sag). 

If Eq. (7) is substituted Into Eq. (6), we obtain an Im- 
plicit relationship between the wavespeed c, the wavelength 
t = 2ir/a, and the flow velocity U. This equation Is written 
out below in dimensionless form: 


1 /-L 4. 2b \ 

¥ (sb * T) 


+ 4r 



2 



C** - U 




F(n) = 0 


(9) 


The following notation has been Introduced 


|D 

r 

D 

D 

r 

r 
I G 


u = pb/m 

r = EhbV(2ir)'* R^D 
c^ S 2 it (D/mb 2 ) *'2 
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c* = c/c 

o 

U* H U/c (11) 

o 

When U* or p Is zero, Eq. (9) has two real solutions for 
the dimensionless wavespeed c*. Above a certain value of U*, 
however, the wavespeeds are complex, and one of the two corresponds 
to unstable motion, that is, waves that grow in amplitude. The 
value of U* at which c* becomes complex depends on the normalized 
wavelength t/2b; the emalleet value of U* is the one that is 
physically significant, since it corresponds to the initial ap- 
pearance of instability. This value of U* is denoted by 
and is called the critical or flutter velocity. The correspond- 
ing values of c* and i/2b are denoted by c*_ and (£/2b) These 
are the critical or flutter wavespeed and wavelength, respectively. 

In Pigs. 3, 6 and 7 of Dowell (1966), c* , U/2b) and U» 

cr cr cr 

are displayed graphically for a flat plate as functions of Mach 
number M and mass ratio p: 


c 


« 

cr 


c* (u,M) 


(V2b)^^ “ ^^/2b)^^ (u,M) 


U* 

cr 


U* 

cr 


(u,M) 


For a flat plate with no pressure differential, Dowell 
showed that the critical wavelength is closely approximated by 
the value which gives the minimum wavespeed in a vacuum (y = 0). 
Setting y « Ap » r * 0 in Eq. (9), we obtain, for the in vaauo 
wavespeed of a flat plate. 
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|c»| 




Inspection of this equation shows that the minimum wavespeed Is 

|c*. I ■ 1 when £/2b ■ 1. As shown In Pig. 6 of Dowell (1966), 
mxifi 

(4/2b) T 1 for any y and M, and especially for M > 1 and 

w 

y << 1. Thus, Dowell*s results for U* (y,M) and c* (y,M) may 

cr cr 

be generalized to apply to a curved plate with a static pressure 
differential by Interpreting them as U* /c , and c* /c ^ , 
respectively. Implicit In this generalization Is the assumption 
that P(n) can be treated as a constant in Eq. (9). Fortunately, 
F(n) does not vary rapidly with the parameters of Interest, i.e., 
y, M, and t/2b. 


As a result of this approximation. Pigs. 3 and 7 of Dowell 

(1966) may be used directly to obtain U* /c . and c* /c . for 

^ cr min cr rain 

a curved plate exposed to a static pressure differential. The 
flutter wavelength is assumed to be the minimum wavespeed in a 
vacuum, and the influence of the pressure differential and the 
cross-stream curvature on the flutter velocity is determined 
entirely by their effect on the minimum wavespeed. 


It remains, therefore, to calculate the minimum wavespeed 
in a vacuum. This Is done by setting y = 0 in Eq. (9), and 
solving for c*: 



(12) 
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The dimensional wavespeed c » c^c* obtained from this 
equation Is shown plotted In Pig. 12, using the following values 
for the dimensions and the mechanical properties of the heat 
shield: 

h - 0.025 cm (0.01 in) 
b » 2.29 cm (0.9 In) 

i^hTin » 305. m/s (1000 ft/sec) 

These values were selected as being representative of the 
several heat shields used in the AEDC tests. They were obtained 
from data supplied by NASA MSPC (see Appendix A). In calculating 
c^ and r, we assumed that the plate bending stiffness D could be 
approximated by its value for a homogeneous elastic plate. 


12(1-vM 

in spite of the heat shields* non-homogeneous construction.* 

This approximation was required because no Independent measure- 
ment of D was available. 

In Pig. 12, the in vacuo wavespeed is shown as a function 
of £/2b for several values of Ap. When Ap * 0, there are two 
values of il/2b that provide the same minimum wavespeed. Of 
these, the longer of the two wavelengths Is found to result in 
the lowest critical velocity U* , because P(n) Is a monotonically 
decreasing function of n» and larger l/2b means emaller n (see 
Eq. [8]) and thus larger P(n) in Eq. (9). By the same reasoning. 


*Some of the heat shields were constructed of nylon mesh Im^ 
pregnated with rubber. 
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FIG. 12. IN~VACUO WAVESPEED IN 3X SCALE MODEL HEAT SHIELD 
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the local minima in the wavespeed c in Pig. 12 for which t/2b >1 i 

remain critical as Ap is Increased. The critical wavelength 
(t/2b)^^ is displayed in Fig. 13 as a function of Ap. Note that I 

increasing Ap causes a rapid deereaae in the critical wavelength, 

y 

and an ineveaBe in the critical wavespeed. ! 

L 

The flutter velocities calculated as described above are 
displayed in Pig. 14. These flutter boundaries are shown va. 

Mach number for several values of Ap. The unstable region in 
the figure is to the right and above these curves. Also shown 
in the same figure is the air flow velocity in the wind tunnel, 
which operated at a constant stagnation temperature of 311®K 
(560°R). The Mach number on the horizontal axis is the local 
Mach number over the face of the heat shield, not the freestream 
or test section Mach number. Recall that the local Mach number 
depends on the angle of attack; as the angle of attack Increases, 
the local Mach number drops from toward 1 (see Pig. 8). 

Thus, at low angle of attack the flow velocity is high; as the 
angle of attack increases, the flow velocity decreases, follow- 
ing the wind tunnel curve in Fig. l4 from right to left in the 
direction of decreasing local Mach number. Eventually, a stability 
boundary is crossed and flutter is auppreaaed. The local Mach 
number at which this occurs is shown in Pig. 15, which is merely 
a cross-plot of Fig. 14. The angle of attack at which the 
stability boundary is crossed will, of course, depend on the free 
stream Mach number in accordan^^e with Fig. 8. 

The AEDC tests were conducted at M_ * 2.75 and 3 .5. There- 
fore. Fig. 14 implies that for Ap ^ 13,790 N/m* (2 psi), flutter 
always occurred at low angles of attack, but was eventually 
suppressed at higher angles. As we shall see, this conclusion 
will be altered when the stability boundaries in Fig. 14 have 
been modified to Incorporate the high Mach number flutter theory. 
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FIG. 14. STABILITY BOUNDARIES FOR 3% SCALE HEAT SHIELDS 
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Critical wavespeeds and frequencies are shown In Figs. 16 
and 17 for Ap * 0 and 13»790 N/n‘ (2 psl). The pressure dlf~ 
ferentlal Ap Increases the critical wavespeed at any given Mach 
nuBber by the sane proportion that It Increases the in vacuo 
wavespeed. The flutter frequency Is given by 


The pressure dlffex*entlal Increases the critical wavespeed and 
at the same time decreases the critical wavelength (see Pig. 12); 
thus, the frequency Is Increased by a greater proportion than 
is the wavespeed. 


A. 3 High Supersonic Theory 

For M > 1.5 or so the assumption of a traveling wave 
solution (Eqs. [A]) is no longer appropriate (Dowell and 
Ventres, 1970). However, a stability boundary based upon the 
use of piston theory aerodynamics and the assumption that the 
length of the plate In the screamwise direction is much greater 
than Its width, can be derived from the theory presented by 
Dowell and Ventres (1970). This stability boundary was derived 
for rectangular plates, so we are retaining the concept of an 
equivalent straight heat shield as discussed previously. 

The stability boundary Is given below. 






(13) 


whe”e B = (M*-l) and q = 1/2 pU* is the local dynamic pressure. 
Note that the stability boundary Is defined in terms of the 
dynamic pressure rether than the velocity as It was at low supc.r> 
sonic speeds. Critical values of q/B calculated from this 
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CRITICAL HAVESPEED VS LOCAL MACH NUMBER 
FOR 3X SCALE HEAT SHIELDS 
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equation for several values of Ap are shown In Table The 
mechanical properties listed on page 30 were used for the com- 
putations. 

TABLE 3. (<*/8)fiutter 


G 

G 

G 




n 

n 

D 





Ap 

(**^®)flutter 

0 

0.139 kN/m* (.0201 psi) 

7 kN/m* (1 psl) 

7.93 kN/m* (1.15 psl) 

13.8 kN/m* (2 psi) 

13.9 kN/m* (3.08 psl) 

27.6 kN/m* (i» psi) 

58.5 kN/m* (8.1»9 psi) 


In Figs. l8 and 19, stability boundaries for = 2.75 and 
3.5 are displayed along with the range of local q/6 encountered 
at each angle of attack. The unstable regime Is toward the top 

of each figure, i.e., > (*1/® Aflutter’ 

tests, the tunnel dynamic pressure was varied at constant by 
changing the tunnel total pressure p^. The dashed curves in 
both figures are the values of q/B that correspond to the maximum 
and minimum tunnel stagnation pressures used at each test section 
Mach number. These dashed curves are terminated at the angle of 
attack where the local Mach number dropped to 1.5. Below this 
Mach number, or alternatively, above the corresponding angle of 
attack, the stability boundaries do not apply, and the low 
supersonic theory must be used instead. 

Both figures show that flutter is suppressed by a pressure 
differential of 28 kN/m* psl) at all tunnel stagnation pres- 
sures, except for a small range of angles near 120° at * 3.5. 
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The stabilizing effect of Ap is much more pronounced at high 
Mach numbers than at the lower speeds. 

In Pigs. 20 and 21 are shown composite curves, for M^ = 2.75 
and M » 3*5, which incorporate both the low and high Mach number 

QO 

stability boundaries. They are plotted as flutter velocity vs. 
local Mach number. Curves are shown for Ap * 0 and Ap * 7 kN/m* 

(1 psi). The stability boundaries for each value of Ap bifurcate 
at M = 1.5; as explained previously, above that Mach number flutter 
occurs at constant q/B rather than constant velocity. The dashed 
portions of each curve are sketched in to provide smooth transi- 
tions between the two Mach number regimes. Also shov;n in both 
figures is the wind tunnel flow velocity vs. local Mach number. 
Flutter is indicated whenever the wind tunnel curve lies above 
the appropriate stability boundary. 

Each figure is for a specific test section Mach number; the 
local Mach number over the face of the heat shield varies as the 
angle of attack changes. Note that at both M^ = 2.75 and 
M^ * 3*5, a pressure differential of 7 kN/m^ (1 psi) suppresses 
flutter over all but a very limited range of local Mach numbers 
when the tunnel is operated at reduced total pressure. Larger 
Ap's would have the same effect at higher tunnel total pressures. 
Since local Mach number is related to angle of attack (recall 
Pig. 8), flutter occurs over only a limited range of angle of 
attack as well. At both M^ = 2.75 and 3*5, flutter occurs when 
the local Mach number is between 1.3 and 2. From Fig. 8, the 
corresponding ranges of angle of attack are 107° < a < 118° at 
M » 2.75» and 116° < a < 126° at M *3-5. Both of these ranges 
agree fairly well with the angles over which the highest fluctuat- 
ing pressure levels were recorded in the AEBC tests. The noise 
spectra measured at these angles show a pronounced low frequency 
content (below 500 Hz) that can reasonably be attributed to the 
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FIG. 21. COMPOSITE STABILITY BOUNDARIES FOR 3% SCALE HEAT 
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heat shield oscillation frequency and its overtones. In this 
respect, therefore, the theory correlates well with the experi- 
ment . 

! The critical Mach number range, where flutter is most likely 

I to occur, lies at the transition between the low and high super- 

I sonic regimes (the dashed portion of the flutter boundaries in 

I Figs. 20 and 21). A single flutter theory, valid for all super- 

i sonic Mach numbers, would be therefore highly desirable. Such 

I j a theory was not attempted, however, because of the numerical 

i L . complications involved. The simplified low and high Mach number 

I theories presented in Sections ^.2 and ^.3 are adequate to explain 

' ‘ the wind tunnel data and to provide guidance for improving the 

I design of the full scale SRB heat shield. 
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5. FULL SCALE SRB HEAT SHIELD: REQUIREMENTS FOR FLUTTER PRE- 

VENTION AND PROPOSED DESIGN MODIFICATIONS 

5.1 Present Heat Shield Design 

The SRB heat shield consists of 2^* segments laced together 
along the radial edges to form a complete annulus. The heat 
shield is attached to the nozzle on the inner perimeter and to 
the aft skirt on the outer perimeter. The segments are sized 
so that each assumes a sag radius of 37.3 cm (14.7 in). 

The segments are made up of three layers of material, sewn 
together in a square quilt pattern. The outermost layer is 
heat resistant #581 Astroquartz woven cloth, and the inner layer 
is #1582 Fiberglass. Between these is sandwiched a layer of 
Piberfrax Insulation. Most of the strength of the heat shield 
derives from the inner and outer layers. 

The surface mass density of the heat shield is about 1.88 
kg/m^ ( 0.012 slugs/ft*). All other mechanical properties of the 
heat shield bearing on its aeroelastlc stability are unknown. 

For most purposes, the heat shield is a limp membrane-like 
structure. With regard to flutter however, its bending stiff- 
ness, however small, is of crucial importance.* The bending 
stiffness of a homogeneous elastic plate is related to the 
Young's Modulus and to the thickness h by the familiar formula 
D = Eh*/12(l-v* ) . The heat shield, however, is built up of 
three different materials and because of its quilted construc- 
tion, the thickness h is not constant. Thus, D must be deter- 
mined by experiment. 


*The bending stiffness D appears implicitly in Eq. (9), which 
determines the flutter velocity at low supersonic Mach numbers, 
and explicitly in Eq. (13) » which determines the flutter dynamic 
pressure at hign supersonic Mach numbers. 
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Also significant with regard to flutter is the ratio be- 
tween In-plane tension and in-plane strain. For a homogeneous 
plate this ratio is Eh. For the fabric heat shield, whereas the 
product Eh has the same physical significance, E and h are not 
separately defined. Eh also must be determined by experiment. 

Because the bending stiffness of the full scale heat shield 
is unknown, we were not able to calculate flutter boundaries 
analogous to those presented in Sec. 3 for the model heat shields. 
Instead we have calculated, using the same theory, the bending 
Btiffneae required to prevent flutter during reentry. The pur- 
pose here is to establish whether or not the heat shield as 
presently conceived may reasonably be expected to have the 
required stiffness. The conclusion reached is that it cannot. 


5.1.1 Stiffness Required to Pievent Flutter at 
Low Supersonic Speeds 

In Fig. 22 is shown the minimum »vavespeed in vacuo for the 
full scale heat shield. The parameter in the figure is the in- 
plane stiffness Eh. These curves were calculated from Eq. (12), 
using b = 76.2 cm (30 in.) and m = 1.88 kg/m^ (0.012 slugs/ftM. 
The in-plane stiffness increases the wavespeed by virtue of the 
added he p stiffness due to the cross-stream curvature R. The 
Increase in wavespeed is only 50% when Eh Increases from 0 to 
17500 kN/m (10® Ib/in. ). 


The local flow velocity over the heat shield during reentry 
is displayed in Pig. 23 for several values of freestream or re- 
entry Mach number. For all local Mach numbers, the highest 
local velocity is encountered when the reentry Mach number is 
3 . 5 . Assuming highest Mach number at which 

the low Mach number theory is accurate), we obtain U = 6'40 m/sec 
(2100 ft/sec). 
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In Fifn. 2H Is shown the maximum local density encountered 
during reentry, displayed as a function of the local Mach number. 
At ■ 1*5 » Piocai * 0,618 kg/m* (1.2 • 10“* slugs/ft*), so 

the density parameter u for the heat shield Is 



0.012 

* 0.25 


I 


I 


I 1 : 

i 1 . 


: 1 



Prom Pig. 7 of Dowell (1966), using M = 1.5 and \i = 0.25, 
we obtain for the critical flow velocity 

U 

— ^ “ 2.5 
^mln 

solving for and using » 6^10 m/sec as calculated above, 

we get 

c . = 268 m/sec (880 ft/sec) 

min 

Prom Pig, 22, the required stiffness is, assuming Eh = 0, 



D » 1980 N-m (li»60 ft-lb) 

If we assume Eh * 1750C kN/m (10* Ib/ln), we get 
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5.1.2 Stiffness Required to Prevent Flutter at High 
Supersonic Speeds 

For M > 1.5» the stability boundary is given by Eq. 13. 

This equation contains (^/B^local l^^an 

Eq. 13 is plotted as q/B vs. D for Ap « 0 and 28 kN/m* (4 psi). 
Note that, when Ap > 0, one can suppress flut«,er by decreasing 
the bending stiffness D. This is not possible when Ap » 0, how- 
ever, and is not possible for any Ap at loim supersonic Mach 
numbers. Making the heat shield more flexible is not, therefore, 
a feasible flutter "fix." 

Figure 26 displays the maximum value of ^^/^)local 
tered during reentry as a function of reentry Mach number. (At 
all reentry Mach numbers, this maximum occurs at the angle of 
attack at which Mj^pcal * The peak value of (q/3)iocal 

from the figure, 

(q/6)iocai “ kN/m^ (2170 Ib/ft*) 

From Fig. 25 » assuming conservatively that Ap * 0, the 
required bending stiffness is 

D = 226c N-m (1670 ft-lb) 



This required stiffness is greater than either of the values 
calculated for ^ 1.5» so it is the controlling stiffness 

requirement for prevention of flutter during reentry. 


To express this bending stiffness requirement in more easily 
comprehensible terms, we have calculated the thickness of a 
homogeneous aluminum plate having this bending stiffness. As- 
suming for aluminum E ■ 6.9 • 10 ’ kN/m* (lO’ psi), the thick- 
ness is 
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REENTRY MACH NUMBER 


FIG. 26. MAXIMUM LOCAL VALUE OF q/6 DURING REENTRY 
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h 


12(l-v^)D 

E 


V, 


- 0.710 cm (0.28 in) 


This is a substantial thickness. The estimate of stiffness re- 
quired is expected to be somewhat conservative, because in de- 
riving Eq. 13, aerodynamic damping was neglected (i.e., y * 0), 
whereas in fact, u = 0.25. 

Nevertheless, the fabric heat shield cannot be expected to 
be as stiff as an aluminum plate 0.710 cm thick. As a result, 
we predict that the heat shield will flutter during reentry. 


5.2 Concepts for Flutter Suppression 

The 3 % scale model studies (both theory and experiment) sug 
gest that a static pre.'sure differential or other means of 
Introducing In-plane tension into the heat shield is beneficial 
in suppressing flutter. The theory also suggests that increas- 
ing the bending stiffness is helpful at both low and high Mach 
numbers. Because of the requirement that the rocket nozzle be 
swiveled during ascent, the amount of radial stiffness that can 
be allowed is limited. We are led, therefore, to consider 
orthotropic structures with greater circumferential than radial 
bending stiffness. All the flutter suppression concepts that 
have been considered to date include one or both of these two 
ingredients (i.e., in-plane tension and orthotropic bending 
stiffness) in some form or other. Before discussing specific 
designs, we will outline briefly the modifications required in 
the low and high Mach number theories to handle orthotropic 
bending stiffness (both theories already incorporate Ap). 
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All that Is required at low Mach numbers Is to modify Eq. 12 
for to Include orthotropic bending stiffness. The aero- 

dynamic theory remains unchanged. Introducing D and D as the 

X y 

circumferential and radial bending stiffnesses, the square of 
the minimum wavespeed becomes 


'min 


2ir' 






(14) 


By suitable choices of D , D and N , N (which may be 

X y X y 

introduced by a pressure differential or by other means), one 
can make sufficiently large to prevent flutter. An applica- 

tion of this formula to a specific design concept Is discussed 
In Appendix B. 


The high Mach number theory > 1.5) is based on the 

low aspect ratio (large length/width ratio) flat plate theory 

presented in Dowell and Ventres (1970). This can easily be 

extended to include D , D , N and N by changing the definitions 

X y X y 

of certain terras in the equations presented there. 


5.3 Proposed Design Modifications 

Desirn concepts considered so far, as well as others recom- 
mended as deserving additional study are listed below, along 
with comments where appropriate. 


(1) Present Heat Shield Vith In~Plane Tension 

The requirement that the heat shield detach from the aft 
skirt at water Impact presently limits the maximum acceptable 
in-place tension to the equivalent of 2 kN/m^ (0.3 psl). Unless 
this limitation is removed (as, for exampel, by using pyrotechnics 
to actively sever the heat shield prior to splash-down), in-place 
tension, by itself, does not merit further consideration. 
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(2) Ppaeant Heat Shield With In^Plane Tension and Ciroumferential 
Stiffness via Civeumferential Hoops 

Setting D ■ 0 In Eq. (H) (zero radial stiffness), we obtain 



This equation has been Incorporated into a design procedure 
(see Appendix B) for selecting the required stiffness and In- 
plane tensions N^, (expressed as an equivalent Ap). The 
required tension Is apparently larger than desirable from a 
structural standpoint, but on the other hand, the theory used 
to establish the required tension may be quite conservative. 

This concept deserves further consideration. 

(3) Present Heat Shield With Orthotropic Bending Stiffness in 
Both Hadial and Cirawnferential Directions via Circumferential 
Hoops and Radial Stays 

The governing equation for the wavespeed is Eq. (l4). No 
detailed calculations have been made as yet, but it is apparent 
that some tension must be retained in the fabric to prevent 
flutter of the fabric between the stiffening members. 

(4) Honconventional Structures , Such as an "Inner Tube" 

Placed Beneath the Fabric Heat Shield 

The pressure required to stabilize the heat shield will be 
quite high, especially if the heat shield is bulged outward during 
reentry. Thus, the "inner tube" may have to be inflated (or the 
pressure increased) after the SRB separates from the Shuttle 
Orbiter. The mechanical arrangement must be such that tension 
is Induced in the heat shield Itself (not merely in the inner 
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tube) to prevent ’’local flutter" of the heat shield (l.e., flap 
ping against the wall of the Inner tube). 

This concept requires further study. Bulging the heat 
shield outward Introduces an additional unknown aspect to the 
problem, because the local flow field will no longer closely 
resemble that encountered In the wind tunnel tests completed 
so far. 
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6. RECOMMENDED PROGRAM OF EXPERIMENTS 

In this section, we outline a program of vibration tests and wind 
tunnel tests that are required to determine how the heat shield 
can be modified to avoid flutter during reentry. These experi- 
ments are needed because the theory, by itself, does not afford 
an acceptable level of confidence In determining whether or not 
a specific heat shield design will flutter. The theory does, 
however, provide valuable insight into the mechanism of the in- 
stability, and helps to identify promising flutter suppression 
concepts. 

A detailed proposal for wind tunnel aeroelastlc tests at 
6 % scale has already been submitted to ?JASA MSFC.* In this 
section, therefore, only salient features of the proposed program 
will be discussed. 


6.1 Aeroelastlc Scaling Parameters 

A formal analysis of the heat shield for the purpose of 
identifying scaling parameters for aeroelastlc and acoustic tests 
has been carried out, and is reported in Volume II of this report. 
The appropriate scaling parameters also appear, of course, in the 
flutter theory presented in Section The scaling parameters 
are listed below (L denotes a typical length). 


M EhL* 
-TT- » 




N L* 

» D ’ 




D 


m 


The significance of each as regards the design ana < onstruc- 
tlon of an aeroelastlcally scaled model and the conduct cf the 
wind tunnel tests will be discussed in turn. 


*BBN Proposal Number P77-NCT-70, Nov. 1976, revised Feb. 1977, 
as P77-NCT-70a. 
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M is the Mach number of the flow over the face of the heat 
shield. According to the wedge flow analogy discussed In 
Section 2, the local Mach number Is the same at all points on 
the heat shield, and is related to the free stream (reentry or 
wind tunnel test section) Mach number as shown In Fig. 8. How- 
ever, this model of the flow adjacent to the heat shield is 
supported by a limited amount of data, and Is in any event a 
simplification of the situation that actually exists. To model 
the flow field a^'-urately, the wind tunnel tests should be con- 
ducted by mounting the aeroelastically scaled heat shield models 
on a geometrically scaled model of the entire aft end of the SRB. 

EhL^/D is a ratio of two structural parameters of the heat 
shield, the in-plane stiffness Eh and the bending stiffness D. 

This parauneter is not very important in determining the stability 
of the present heat shield design (see Section 5.1.1), but might 
become significant if the design were modified substantially. 

In any event, EhL^/D can be matched between model and full scale 
by proper selection of the material of which the model is con- 
structed. 

N^L^/D and NyL*/D are the in-plane tension loading parameters. 
In-plane tension can be introduced by imposing a pressure differen- 
tial Ap across the heat shield, or by mechanical means. If the 

tension is solely due to Ap then N = ApR/2, and N « ApR, where R 

A y 

is the sag radius. In this situation, the two parameters reduce 
to one. Because R ~ L, we may write this one parameter as ApL*/D. 
In a wind tunnel test, Ap can be set at any desired value by con- 
trolling the pressure within the cavity behind the heat shield. 

If and/or are inti’oduced by mechanical means at full scale, 
similar methods may be used in the model. Tlius, these parameters 
impose no conceptual difficulty in scaling. 

The remaining two scaling parameters, qLVd and pL/m, deter- 
mine the bending stiffness and surface mass density required of 
the model heat shield. Because most wind tunnels operate at a 
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fixed total temperature, the test section flow density u... 
dynamic pressure cannot be separately controlled. Therefore, 
while both of these parameters can be scaled at any single test 
section Mach number by appropriate choices for D and m, different 
values of D and m will be required at other Mach numbers. In con- 
sequence, perfect aeroelastlc scaling Is not possible at all Mach 
numbers. The heat shield models can be constructed, however, so 
as to make the scaling conservative (i.e. so as to decrease the 
dynamic pressure at which the model flutters). Theory indicates 
that the dynamic pressure at which flutter occurs generally in- 
creases as the parameter pL/m Increases. Conservative scaling 
would be achieved, therefore, if ( £ ^^^^"'^full scale 
at all test section Mach numbers. 

6.2 Vibration Tests 

A series of vibration tests on both the scale model heat 
shields and on the full scale heat shield are recommended prior 
to wind tunnel entry. Such tests are recommended in NASA Space 
Vehicle Design Criteria Document on Panel Flutter, NASA SP-8004. 
The purpose of these tests is to verify that the wind tunnel model 
is dynamically similar to the full scale article, and to verify 
that the theoretical flutter model adequately represents the 
in-vacuo vibration behavior of both the model and the full scale 
structure . 

These tests will be carried out by shaking the heat shield 
locally at various frequencies, and measuring the amplitude and 
speed of the waves that propagate around the annulus. If a ccm- 
plete full scale heat shield (24 segments) is not available, 
meaningful results could probably be achieved with only 4 or 5 
segments . 
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6.3 Wind Tunnel Tests 

The wind tunnel tests should be conducted at test section 
Mach numbers between 2 and 3.5> because this Is the most critical 
range with regard to flutter. A range of dynamic pressures should 
be employed, with the scaled maximum f^ynamlc pressure reaching or 
exceeding a level 5055 greater than the maximum levels anticipated 
during reentry (Fig. 2). The igle of attack range should also 
encompass the range expected during reentry (Fig. 3), but tae 
highest priority should be assigned to those angles at '^ac'' test 
section Mach number at which the local Mach number over + ae heat 
shield is between 1 and 2. 

Testing a rigid heat shield, inscrumf^nted ; >r th ? measure- 
ment of local velocities and pressures, over the entire range of 
Mach numbers and angles of attack would clarify the details of 
the flow field, and would facilitate comparisons between theory 
and experiment for the aeroe.lastlc tests to follow. 

It would be desirable, but not absolutely essential, to 
provide means for varying the angle of attack between 90° 
and l80° without shutting down the tunnel. This would allow 
the tunnel to be started at a reduced dynamic pressure with the 
angle of attack set at l80° (heat shield facing the flow), where 
flutter is not anticipated. Two test procedures are then possi- 
ble. The model can be rotated to the desired angle of attack, 
and then the dynamic pressure increased until flutter is encoun- 
tered, or until q is reached. This procedure would assure that 
the flutter boundary is encountered only once. Alternatively, 
the dynamic pressure can first be raised to the desired val'.e, 
and then the model can be rotated through the desired angle of 
attack range, Depending upon the heat shield design, two flutter 
boundaries may be encountered, one where flutter commences, 
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another where flutter ceases. This is because the local Mach 
number passes through the most critical range (1 < < 2) at 

Intermediate values of a Csee Pigs. 8, 20, and 211, 


Measuring the motion of the model heat shields presents some 
difficulty, because the amplitude of the expected oscillations are 
on the order of the width of the annulus, and the heat shields 
will be very light In weight. An optical sensing system that 
adds no mass whatsoever to the heat shield Is described In the 
proposal. High speed films will be useful In studying the 
details of the oscillations, as they were In the AEDC tests. Other 
Instrumentation for measuring flow velocity and dynamic pressure 
will be conventional. 


Several flutter suppression concepts should be tested, 
although not all need be tested over the entire range of M, q, 
and o. Purtheremore, three versions of each "fix" should be 
constructed and tested. These will have graduated levels of struc- 
tural stiffness (or some other key parameter, such as Internal 
pressure in the case of the Inflated "inner tube" concept). This 
procedure will mak it possible to interpolate or extrapolate the 
data to determine the stiffness level required to suppress flutter 
during reentry. The three levels of stiffness need not correspond 
to values that would be acceptable in full scale; it would be more 
informative to discover the actual stiffness required to prevent 
flutter, even though this level be unacceptable on other grounds, 
than to complete the test program without having tested a flutter- 
free heat shield. 
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7. SUMMARY AND CONCLUSIONS 

The oscillations observed in the AEDC wind tunnel tects on 
a 3 % scale model of the SFu) have been traced to aeroelastic in- 
stability of the heat shield. A flutter theory based on a sim- 
plified geometrical model of the annular heat shield shows 
relatively good agreement with experiment as to the frequency, 
wavelength and wavespeed of the oscillations, and the range of 
angles of attack at which the oscillations are most severe. 

It was not possible to apply the same theory directly to the 
full scale SRB heat shield to determine whether or not it would 
flutter during reentry, because the requisite mechanical properties 
of the heat shield are not known. Instead, the bending stiffness 
required to prevent flutter was calculated. This stiffness was 
found to be far greater than the fabric heat shield could possib- 
ly possess in its present form. On this basis, the SRB heat 
shield is expected to flutter during reentry. 

Theory indicates that flutter can be suppressed by providing 
sufficient in-plane tension in the heat shield, or by increasing 
its bending stiffness (the present heat shield is a limp fabric 
structure). In-plane tension can be supplied by placing a static 
pressure differential between the inner and outer surfaces of the 
heat shield, or by any equivalent means. The amount of tension, 
acting alone, that would be required to stabilize the heat shield 
would be quite large. Combined tension and bending stiffness is 
suggested as a better means of suppressing flutter. 



A program of experiments, including vibration tests on both 
the full scale heat shield (or a portion thereof), and on aero- 
elastlcally scaled models of the heat shield, and a series cf 
aeroelastic wind tunnel tests is proposed in order to determine a 
design modification or "fix" that will stabilize the heat shield. 
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The wind tunnel tests should be conducted at test section Mach 
numbers between 2 and 3*5. Dynamic pressure should be Increased 
until flutter Is encountered, or until maximum tunnel aynamlc 
pressure Is reached. Maximum tunnel dynamic pressure should 
correspond at model scale to a level 50 % greater than the maximum 
pressure expected during reentry. The tests should be carried out 
on a geometrically scaled model of the aft end of the SRB. This 
model should be rotated through the full range of angles of 
attack expected during reentry, although highest priority should 
be assigned those angles of attack that produce local Mach num- 
bers, over the face of the heat shield, between 1 and 2. This 
Mach number range Is most critical with regard to flutter. Each 
flutter suppression concept should be tested In three different 
versions, having different levels of stiffness. This procedure 
will make It possible to Interpolate or extrapolate the data to 
determine the minimum stiffness required to prevent flutter during 
reentry . 
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APPENDIX A 

ELASTIC CONSTANTS OF THE 3X SCALE MODEL HEAT SHIELDS 
USED IN THE HIND TUNNEL TESTS AT AEDC 

Using data provided by R.W. Walker of NASA MSFC (private 
communication), the In-plane stiffness Eh, and compresslonal 
wave speed «^Eh/m were calculated for the model heat shields used 
In the wind tunnel tests at AEDC. Walker performed uniaxial pull 
tests on rectangular samples of material cut from the heat shields 
at two orientations, one radial, and the other circumferential. 

He supplied graphic recordings of force vs elongation. These 
curves were not linear, so Judgement had to be exercised In 
deciding at what strain level (or elongation) to measure the 
slope. In consequence, the valves calculated for Eh and /Eh/m 
are at best crude estimates of the equivalent linear elastic 
constants of the heat shields. 

1. Materials Used at = 2.75 

A. Lightweight Nylon Impregnated with Rubber: 

Sample Orientation 



Radial 

Circumferential 

h 

.015 

cm 

m 

0.20 

kg/m* 

Eh 

1.8 • 10“ N/m 

2.3 • 10“ N/m 

Eh/m 

300 m/s 

335 ra/s 
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B. Polyester Cloth 

Sample Orientation 




Radial 

Circumferential 


h 


.025 cm 


m 


.097 kg/m* 


Eh 

1.1 • lO** 

N/m 3.7 • 10 * N/m 


Eh/m 

335 m/s 

195 m/s 

Materials Used at » 

3.50 


A. 

Heavy Nylon Impregnated with Rubber: 



Sample Orientation 



Radial 

Circumferential 


h 


.025 cm 


m 


.028 kg/m* 


Eh 

1.5 • lO” 

N/m 1.5 • lO** N/m 


Eh/m 

231 m/s 

231 m/s 

B. 

Polyester Cloth 

Sample Orientation 



Radial 

Circumferential 


h 


.025 cm 


m 


.097 kg/m* 


Eh 

9.3 • 10 ® 

N/m 3.5 • 10 ® N/m 


Eh/m 

309 m/s 

190 m/s 
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APPENDIX B 

HEAT SHIELD STIFFENED BY CIRCUMFERENTIAL HOOPS AND 

IN-PLANE TENSION 


As presently conceived, the heat shield is a limp fabric 
curtain. The curtain covers an annular opening having an 
outer diameter of 208 inches, and an inner diameter of l^l8 
Inches. Enough slack is provided for so that the fabric can 
assume a sag radius of about 15 inches. 

One way to Improve the aeroelastlc stability of the heat 
shield is to stiffen it by 

(1) sewing circumferential rings or hoops into the 
fabric at equally spaced radial increments, and 

(2) applying an in-plane tension to the fabric. 


At low supersonic Mach numbers (1 < M < 1.5), the criteria for 
stability is that the minimum in-vaouo wavespeed exceed 268 m/s 
(880 ft/sec). (See Section 5.1*1.) The minimum wavespeed in the 
reinforced heat shield is given by Eq. (15) , 



and the wavelength at which this minimum wavespeed occurs is 



If the inplane tensions and N are caused by static pressure 
differential, then N = N„/2 (see Eq. [3]). Therefore, 

X J 



JL 

2m 


68 



I ^ I 


Bolt Beranek and Newman Inc. 


expression as follows: I 

1 

1 

. 1/2 ly (B-3, I 

I 

j 

4 

If the critical wavelength greatly exceeds the width of the heat 1 

shield, then the "straightened out" geometrical model upon which f 

these calculations are based (Section 4.2) becomes Inconsistent 
Hence restrict t/2b to be less than, say, 2.11. From Eq. (B-2), this 
Is equivalent to requiring that 



= 2 


O * - l. 
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It Is convenient to re-arrange this 


'min 


2v 


(t?) 





Inspection of Eq. (B-3) shows that large values of D /N b are 

X y 

advantageous. In that they reduce the In-plane tension needed to 
meet the required Thus set 



Eq. (B-3) becomes 


2 . 


°mln ‘ (j) 


(B-4) 


Solving for the required 


N 

_2.» 

m 
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!!*■ 

m 


°mln 

2¥7T^ + 1/2 


(B-5) 


Prom Eq. the associated circumferential bending stiffness 

Is 



(B-6) 


The circumferential bending stiffness Is due entirely to 
the hoops. If N hoops are efnployed, each having the same beam 
bending stiffness El, then 


On the other hand, the mass/unit area m depends on both the 
surface mass density of the fabric and the mass of the hoops: 

” ’ "fabric 

s 

where is the surface area of the heat shield, taking the sag 
radius into account. 

If we use tubular hoops with wall thickness t, then 
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I ■ wr t 


NEitr*t 


X b 


Solving for r. 


(V ] 

y NEirt f 


Let L be the average circumference of the hoops and p. be the mass 

h 

density of the material of which they are constructed. Then their 
total mass Is 


M. * 2. TtLp^ 

hoops h 


= 2Nit ( ^ \ ‘''3 tLp 
\ NExt / ” 


So the mass/unit area of the reinforced heat shield is 


m * m. ^ . +2Nir 
fabric 


(V ) 

XNETTt / 


tLp^/Ag 


(B-7) 


Eqs. (B-6) and (B-7) can be solved Iteratively for the mass/unit 

area m and the required bending stiffness D^. The procedure is 

to assume a value for m (m * m„ ^ ^ is a good starting assump- 

fabrlc 

tion), compute D from Eq. (B-6), then compute a new value for 
m from Eq. (B-7), and so on. 


For example, if we use four aluminum tubes with a wall 


thickness of 0.127 cm (0.05 in), and assume that m 
kg/m^ (.012 slugs/ft*), then 


fabric 


1.88 
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N - it 

7.10’ kN/m\ do’ psi) 

277 kg/m* (.00311 slugs/in*) 


and we obtain 


m = it. Oil kg/m’ (0.0257 slugs/ft’) 
D » 3.59-10‘‘ N-m (2.65*10- ft*lb) 

X 


The radius and wall thickness of the it tubular aluminum hoops are: 


r = 2.9 cm 
t “ 0.13 cm 


(1 15 in) 
(0.05 in) 


The required In-plane tension N is, from Eq. (B-5), 

N * 20.0 kN/m (2120 Ib/ft) 

y 

The pressure differential required to establish this tension is 

Ap = jr 
R 


= 8l.it kN/m’ (11.8 psi) 
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K 

8 

b 

c 

°mln 

Co 

D 

E 

Eh 

f 

P(n) 

h 

1 

M 
M 


local 


M 


m 


N 

N 

P 

q 

R 


X 

y 


LIST OF SYMBOLS 

“ surface area of heat shield 
« Width of annular heat shield (Pig. 11). 

■ Wavespeed. 

■ Minimum wavespeed. 

■ Minimum in~vaeuo wavespeed In a flat plate. 

= Bending stiffness. 

® Young's modulus. 

“ In-plane stiffness. 

= Frequency. 

» Aerodynamic Integral In Eq. (9) (from Dowell [196^]). 
= Thickness of model heat shield. 

= 

= Wavelength. 

= Mach number. 

* Mach number over face of heat shield. 

= Test section or reentry Mach number. 

■ mass/unit area of heat shield. 

» Stream-wise in-plane tension. 

« Cross-stream In-plane tension. 

“ Pressure. 

*s Dynamic pressure (^ pU*) 

■ Sag radius of heat shield (Pig. 11) 
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t ■ Time . 

V “ Plow velocity. 

w ■ Transverse deflection of heat shield. 

X » Stream-wise coordinate (Fig. 11). 

y » Cross-stream coordinate (Fig. 11). 


a 

6 

r 

Ap 

n 

y 

V 

p 


- Wavenumber (2it/]l) 

= / M^-1 

* Dimensionless constant defined in Eq. (10). 

= Static pressure differential. 

= Dimensionless variable defined in Eq. (8). 

= Normalized flow density (pb/m). 

= Poisson’s ratio. €, 

= Flow density. 

* Airy stress function in Eq. (2). 


Superscri pts 

* 

* Flow velocity or wavespeed d vided by c ^ (minimum 
in-vacuo wavespeed in a flat plate). 

Subscripts 

cr ■ Critical value or value at which flutter occurs, 
min * Minimum value. 

«> ■ Free-stream value . 


w • > 


7^4 
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